Semiclassical dynamics of domain walls 
in the one-dimensional Ising ferromagnet in a transverse field 
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We investigate analytically and numerically the dynamics of domain walls in a spin chain with 
ferromagnetic Ising interaction and subject to an external magnetic field perpendicular to the easy 
magnetization axis (transverse field Ising model). The analytical results obtained within the con- 
tinuum approximation and numerical simulations performed for discrete classical model are used to 
analyze the quantum properties of domain walls using the semiclassical approximation. We show 
that the domain wall spectrum shows a band structure consisting of 25 non-intersecting zones. 
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I. INTRODUCTION 

Nonlinear topologically nontrivial excitations, or soli- 
tons, are known to exist in lower-dimensional magnetic 
systems. There is both theoretical and experimental ev- 
idence of solitons and in some cases soliton effects dom- 
inate the thermodynamic behavior of one-dimensional 
(ID) magnets (spin chains). Soliton effects possibly 
emerged first when absence of long range order in lower- 
dimensional magnets was proved: From simple entropy 
arguments it has been shown that long range mag- 
netic order is not possible in one-dimensional magnets at 
nonzero temperature and it became clear later that kink- 
type solitons, which in fact are magnetic domain walls 
(DW), must be considered as elementary excitations at 
nonzero temperature in one-dimensional magnets, for re- 
views seei*£. Whereas solitons in one-dimensional mag- 
nets have not been directly observed, dynamic soliton 
effects such as soliton motion and the soliton-magnon in- 
teraction result in soliton contributions to the dynamic 
response functions, which can be studied experimentally. 
For example, solitons contribute to the specific heat and 
to the linewidth of electron spin resonance and the trans- 
lational motion of kinks leads to the so-called soliton cen- 
tral peak, which can be detected through neutron scat- 
tering experiments 

In most approaches solitons have been considered using 
classical continuum models, such as the Landau-Lifshitz 
equations or the sine Gordon equation, see Refs. Ill2j |il 
On the other hand, in order to describe a material such 
as CSC0CI3 the XXZ-model with spin 5=1/2 has to 
be used, and solitons occur as quantum objectsAa Even 
for the material CSMF3, the well-known standard exam- 
ple for classical solitons in ID magnets, quantum effects 
due to its spin 5 = 1 arc essential. 1 Although a detailed 
analysis concentrated on classical models, considerable 
achievements in the field of quantum solitons were ob- 
tained rather early: We note, first of all, that the first 
nonlinear excitations (spin complexes) were investigated 
by Bethe in the isotropic one-dimensional ferromagnet 



as long ago as 1931, essentially in parallel with the pre- 
diction of magnons£ Currently it is established that the 
more general XYZ-model with spin 5 = 1/2 is exactly 
integrable, and the quantum nonlinear excitations (spin 
complexes) are known from the solution of this modeli 

At first sight quantum and classical solitons constitute 
essentially different objects. The main property of classi- 
cal solitons (or more precisely: solitary waves) is localiza- 
tion, while a quantum soliton is characterized by a defi- 
nite value of the quasi-momentum P and in virtue of this 
it is spatially delocalized. However, this contradiction is 
removed if one investigates the spin deviation localization 
in the coordinate system with the origin moving with the 
group velocity of the solitoni 3 . The comparison of classi- 
cal and quantum solitons^ reveals a striking feature of 
the XYZ-model: the dispersion law (the dependence of 
energy E on the momentum P) of a spin complex in this 
model with spin 5=1/2 exactly coincides with the cor- 
responding dependence found for the soliton in the clas- 
sical Landau-Lifshitz equations. It is evident that such 
an exact correspondence cannot be a general rule; prob- 
ably it is associated with the exact integrability of both 
models. On the other hand, it has become clear that if 
one compares the characteristics which are relevant for 
both quantum and classical approaches, first of all the 
dispersion law E — E(P), then there are no fundamental 
distinctions between quantum and classical solitons. Re- 
newed interest in the problem of quantum properties of 
domain walls was stimulated due to the study of quan- 
tum tunneling DW chirality effects jgaLiiLiiiiSi and also by 
the prediction of new effects of destructive interference in 
kink tunneling between neighboring crystal lattice sitesia 
and Bloch oscillations of the solitons. 13 

The structure and the properties of solitons in the one- 
dimensional Ising model in the presence of a transverse 
magnetic field have been studied some time agoii It has 
been demonstrated in the semiclassical limit of large spin 
values 5 ^> 1 how classical localization correlates with 
quantum derealization. Within the framework of the 
quantum approach the emergence of a band structure 
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has been revealed. However, the quantum kink disper- 
sion law found in this work has not been compared with 
the classical one. In Ref. ^] only static solutions of the 
classical problem have been used, and the comparison 
of the quantum kink dispersion law with those obtained 
within the semiclassical quantization model has not been 
made. 

In this article we investigate analytically and numeri- 
cally the dynamics of domain walls in a spin chain with 
ferromagnetic Ising interaction and subject to a mag- 
netic field perpendicular to the easy magnetization axis 
(transverse field Ising model). Analytical results ob- 
tained within the continuum approximation and numeri- 
cal simulations performed for the discrete classical model 
are used for the analysis of the semiclassical dynamics of 
domain walls including an account of the lattice pinning 
potential. We show that the spectrum of the DW in the 
classical continuum approximation is characterized by a 
periodic dependence of the kink energy on its momen- 
tum. This produces a number of non-trivial features for 
the DW motion in the presence of the pinning potential. 
The quantum properties of the domain wall are discussed 
on the basis of the semiclassical approximation. For the 
semiclassical dynamics the role of lattice pinning effects 
in the formation of the band structure of the domain wall 
spectrum is investigated. The presence of the periodic 
character of the DW dispersion law with lattice pinning 
taken into account results in a band structure for E(P) 
with 25* non-intersecting branches. 

The outline of this article is as follows. In Sec. ITT1 we 
formulate the discrete classical model and study the DW 
properties within this model. The main results of this 
section are the introduction of the DW coordinate, X, 
and the calculation of the pinning potential for the DW, 
U(X). This is then used for the analysis of both the clas- 
sical dynamics of the DW and its semiclassical quantiza- 
tion. In Scc. lIIII the complementary continuum approach 
is used to investigate the motion of the DW with finite 
velocity. The DW linear momentum P and the disper- 
sion law of the DW, E = E(P), are calculated, and E{P) 
is found to be a periodic function of the momentum. The 
results for U(X) and E(P) are employed in the next two 
sections to describe the dynamics of the DW in the frame- 
work of the method of collective variables, assuming that 
the DW is an effective quasiparticle with kinetic energy 
E(P), moving in the potential U(X): In Sec. HVI the 
specifics of forced DW motion as well as the detailed fea- 
tures of the band spectrum of the DW in the case of the 
weak potential U(X) are discussed. The peculiarities of 
the dynamics of the DW in a finite potential associated 
with the periodic character of E(P) are considered in 
Sec. [V] Here it becomes clear that finite motion is typ- 
ical for states which are either close to the minimum or 
close to the maximum of the potential U(X). In Sec. I VII 
we consider the quantum tunneling transitions between 
these states, corresponding to an adjacent unit cell. In 
this section the general character of the DW dispersion 
law including the effects of lattice pinning is discussed. 



Sec. IVIll gives our concluding remarks. 

II. THE DISCRETE MODEL AND A DOMAIN 
WALL STRUCTURE. 

We start from the Hamiltonian describing a spin chain 
with Ising-type exchange interaction in a magnetic field 
H directed perpendicularly to the easy axis of a magnet. 
Spins S n are located at points na of a chain (distance a, 
n integer), 

H=-jJ2S*S* n+1 -gt, B Hj2SZ- (1) 

n n 

J is the exchange integral, g the gyromagnetic ratio and 
Hb the Bohr magneton. The dimensionless field h — 
gfisH/^JS will be used in the following (the notation 
2h = 7 was used in RefiAi) This model is usually called 
transverse field Ising model. 

For classical spins the ground state of the model at 
h < 1 {g^sH < 2JS) is doubly degenerated with 

S x = Ssintfo, S z = ±Scos6» , smO Q = h. (2) 

We begin to describe the semiclassical motion of the 
DW within the discrete spin model by studying the solu- 
tions for moving DW's in the discrete classical model. We 
notice that the exchange part for the model of Eq. (JJJ is 
the simplest particular case of the so-called XYZ-modcl 
with three independent parameters Ji =/= J2 J3 

Hxyz = —^2(JiS£S% +1 + J2S%Sn +1 + J 3 S*S* +1 ). 

For an isotropic FM we have J\ = J2 = J3 and DW's 
strictly are absent (however, 'pulse solitons' may exist) 
For the XXZ-model (Ji — J2 < J3) in the absence of 
the magnetic field the exact static solution for the dis- 
crete model with classical spins has been constructed by 
Gochev^. As far as we know this discrete problem has 
not yet been solved for H 7^ or for the more general clas- 
sical XYZ-model even for H = 0. The solution found by 
Gochev for the XXZ-model reads S* = S t&nh{n(n — no)], 
with k — ln[(J3 + y/ Jf — J\)/Ji\, where X/k measures 
the domain wall thickness. In the Ising limit the DW 
thickness goes to zero, l/ft^OatJi^O. We empha- 
size that the value of n which describes the DW position 
in this solution (taken as n — no above) is an arbitrary 
(not necessarily integer) number. Thus the model with 
H = has the nontrivial property that the energy of the 
DW does not depend on its center position for arbitrarily 
large anisotropy. It is not clear, whether this property is 
valid for the same Hamiltonian in the quantum case. 

For the Ising model we easily understand the above re- 
sult since at H = the DW is described by the following 
solution: S* = —S at n < no, = +S at n > no, how- 
ever *S~| at n = n may have arbitrary values. If SjL = 
the DW is localized on the spin with n = no', on the 
other hand, for S* — —S or — +S at n — n the DW 
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localized in the center of the bond which connects the 
spin at n — n with spins at n = Uq + 1 or n = tiq — 1, 
correspondingly. For H — the energies of DW's cen- 
tered on the spin (central spin (CS) DW) or on the bond 
(central bond (CB) DW) coincide. For all intermediate 
cases S* ^ 0, ±1 at n — uq it is natural to postulate 
that the DW is localized on a point X, which does not 
coincide with a site or the bond center and to describe 
the DW dynamics in terms of its coordinate X treated 
as a collective variable. Let us introduce X in the fol- 
lowing way: we choose some lattice site and define the 
DW located on this site to have the coordinate X = 0. 
Let us then find the z-projection S^J^ of the total spin 
Sf ot of the DW localized on that site. It is then natural 
to determine the coordinate of any DW via the value of 
the total z-projection of the spin S^ ot , associated with 

this DW from the expression X = a(S^ ot — Sl ^)/2S. 
In order to account appropriately for possible divergen- 



cies in infinite chains the difference S% ot — c 



<S tot w is calcu- 



lated as E^KPO " S z n (X = 0)] (here and 
S*(X = 0) are the z— projection of spins on site n for 
the DW with the coordinate X and the DW with X = 
resp.). Then, for the DW in the model of Eq. Q with 
H 7^ the complete definition reads 



^CB'^CS 



Ising 

Single-ion, 1/20 




X = 



a 



n— oo 

K(X)-S* n (X = 0)}. (3) 
2b cos #0 ' — ' 

n— — oo 



For zero field the DW energy E does not depend on the 
coordinate X. At H ^ 0, as shown in Ref. IJ, a de- 
pendence of the DW energy on its coordinate emerges, 
pointing to the presence of a lattice pinning potential 
U = U(X). The pinning potential U(X) has its mini- 
mum value for X = na, i.e. for the CS DW, and its max- 
imum value for X = (a/2)(2n+ 1), i.e. for the CB DWM. 
In order to investigate the nonstationary dynamics of the 
DW we will treat U(X) as a potential energy. The anal- 
ysis of the related kinetic energy will be discussed in the 
next section. 

When the DW is moving in a discrete chain it is natural 
to use the periodic potential U(X), U(X + a) = U(X). 
The form of the potential can be chosen as inii, 



U{x) = U sm 2 (irX/a), 



(4) 



where Uq characterizes the intensity of pinning caused by 
discreteness. With the choice of Eq. (J1J for the potential, 
the value of Uq can be found as the energy difference be- 
tween the static central spin and central bond DW'sii^ In 
order to calculate Ecs and Eqb wg need in the solution 
of the corresponding classical discrete problem, which is 
known for small magnetic field only^ Moreover, for our 
purpose we need to know not only the energy difference 
Ecs — Ecb, but the full dependence U{X) in order to 
verify the dependence of Eq. J2J| . 

To reveal U (X) we have carried out a numerical analy- 
sis of the model of Eq. in accordance with the defini- 
tion Eq. f° r the coordinate X. We have searched the 



FIG. 1: The difference of the energies (in units of JS ) for CB 
and CS domain walls (shown schematically at the bottom of 
the figure) vs. magnetic field. The corresponding dependence 
for the model with single-ion anisotropy, normalized by 20, is 
also given for comparison (open circles). 



conditional minimum of the energy as given by Eq. Q 
for a finite spin chain at a fixed value of S^ ot . To solve 
this problem, the simplex type method with non-linear 
constraints was chosen as the method of minimization. 
This method is based on the steepest descent routine ap- 
plied to functions of a large number of variables, and it 
is able to find the conditional minimum for a given func- 
tion with fixing a small number of combinations to given 
values. The method exhibits a fast convergence when a 
spin distribution with only one domain wall placed near 
the point of inflection is used as starting condition. For 
our problem, the angular coordinates of each spin were 
chosen as variables, and the minimum of the energy was 
found with fixing the value of 5 t z ot . The DW was cre- 
ated by fixing the direction of two spins on the ends of 
the chain corresponding to Eq. @. The chain length N 
varied from 20 to 100, and for the the case of interest, 
H < JS, the result was independent of the chain length 
for N > 30. The value of Uo in Eq. @ was determined 
from the difference Ecs — Ecb', the behaviour of this 
quantity is shown in Fig. ^ 

The analysis shown in Fig. [21 demonstrates that in the 
case of interest to us, U(X) is fairly well described 

by Eq. Q. If we consider a more general dependence, 



4 



/ \ 0.02- 


U(X) 

it V 


V0.01 - 


'/ \ 
'/ V 

'/ V 
7 V 

/ V 



2 XI a 



FIG. 2: Shape of the domain wall pinning potential U (X) (in 
units of JS 2 ) for the dimensionless magnetic field h = 0.05. 
Symbols denote the numerical data. The empirical depen- 
dence 0.024 sin 2 (wx/2a), and a fit for a more general depen- 
dence with Uo = 0.024, Ui = —0.00255 are shown as dashed 
and full lines, respectively. 



allowing for one more Fourier component, 

U{x) = U sin 2 (7rX/a) + U x sin 2 (2vrX/a), (5) 

the dependence of U(X) is reproduced with a deviation 
of less than 0.1%. However, as the correction related 
to U\ is small, we will use in the following mainly the 
simplest expression, Eq. (0J. 

Note, as Uo — » for h — *■ 0, and for h — *■ 1 the trans- 
verse field Ising model is characterized by a rather small 
pinning potential. For comparison, in Fig. ^ the cor- 
responding dependence is presented for a ferromagnet 
with an isotropic Heisenberg interaction and a single ion 
anisotropy of the form W a — — ^Z n K{Sf l ) 2 1 which has 
the same amplitude as Eq. 10} at K = J. For this model 
Uq has a maximum at H — 0; the typical values of Uo 
are approximately 20 times higher than for the transverse 
field Ising model. Furthermore, the quantity Uo has op- 
posite sign, i.e., in contrast to the transverse field Ising 
model one has Ecb < Ecs- 



III. DOMAIN WALL DYNAMICS IN 
CONTINUUM APPROXIMATION. 



Most of the results for the dynamics of domain walls or 
kink-type solitons in magnets have been obtained in the 
continuum approach, replacing S n (£) by the smooth vari- 
able spin density S(x, t), where x is the coordinate along 
the chain. Within this approach the dynamics of the vec- 
tor field S(x,t) is described by the well known Landau- 
Lifshitz equation without dissipation^ (see alsoS^) 



SW{S} \ 




Here W{S(x,t)} is the ferromagnetic energy as func- 
tional of the spin density. For our model, the transverse 
field Ising chain, W{S(x,t)} corresponds to the Hamil- 
tonian of Eq. Q and is written as*i^ 



JS Z - gfi B HS x 



Using the continuum approach allows us to find a solu- 
tion, which describes a DW moving with a given velocity 
V. On the other hand, the discreteness effects are evi- 
dently lost going from Eq. Q to Eq. (0). 

We use the relation S 2 = S 2 = const to write S = 
S m(x, t) i.e. to express the spin through the unit vector 
field m(x, t). Its direction is determined by two indepen- 
dent variables, we will use the angular variables 9 and 
f, 



m - 



cos#, 



m, 



sin 9 cos if, m y =sin#sin< / 9 



In terms of these variables the Landau-Lifshitz Eq. © 
reads 



ns . 

— sin 

a 


89 
dt 


SW 

Sip ' 


hs . , 

— sine 
a 


} dtp 


SW 


dt 


~~S9 



(G) 



(8) 



where VFj^, <p} is the ferromagnetic energy written as a 
functional of 9 and tp. These equations can be consid- 
ered as classical Hamiltonian equations for the canoni- 
cally conjugate variables cos 9 (momentum) and p (coor- 
dinate), with Hamilton function W. 

Kink solitons in such a continuum approach have 
been described for a number of magnetic chain models. 
These include the Heisenberg chain with two single ion 
anisotropiesiZiiSiiii, the Ising chain with transverse ex- 
change breaking the xy— symmetrySS and the xy— like 
Heisenberg chain with an external symmetry breaking 
field^li2£ The qualitative results for kink solitons with 
permanent shape are analogous in all these models: The 
dispersion relation, energy vs. velocity, has two branches 
and one distinguishes between the lower energy domain 
wall (LDW) and the upper energy domain wall (UDW). 
The energies of these two solutions merge at the max- 
imum possible velocity V = V c . A more natural for- 
mulation results when the momentum is introduced in- 
stead of the velocity: then the dispersion relation be- 
comes single valued and periodic with the magnetic unit 
cell^ In particular the dispersion relations can be given 
explicitly for the Heisenberg chain with two single ion 
anisotropiesiiii^ 

Eqs. © or © are usually considered as purely classical 
equations. To discuss their applicability to the quantum 
regime one may use the quantum-field approach based 
on the spin coherent states formalism;2A Then the quan- 
tum phase (Berry phase) appears, leading to qualitative 
effects such as the suppression of quantum fluctuations 
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for antiferromagnetic chains 2 ^ and small magnetic parti- 
cles with lialf-odd-integer total spiniSi In this approach 
the spin state on every site n is determined by the spin 
coherent state |m), for which S \m) = Srfi \m). Here m, 
as before, is a unit vector, rh 2 = 1. In this approach, the 
dynamics of the mean value of spin S = Sm is described 
by a Lagrangian which can be written in the form 



£{S} = hS I —A§-W{S}, 



where 



MS) 



x S 



S(Srt + S) 



(9) 



(10) 



ft is a unit vector with arbitrary direction, denoting the 
quantization axis for coherent states. It is important that 
A has the form of the vector potential of a magnetic 
monopole field in the full space {S} (not subject to the 
constraint S 2 — const). 

The vector potential has a singularity (Dirac string) 
for S n = —5, i.e., on a half-line in space {S}. Usually, 
the "north pole" gauge with n = e z is used, then the 

quantity A(S) S acquires the familiar form 



AS = (1- cos 6) {d<p/dt). 



(11) 



In the saddle point approximation for the Lagrangians 
Eq. @ or (|llfl one recovers the classical Landau-Lifshitz 
equations in the form of Eq. © or (JSJ for the mean value 
of spin S. The potential A of the monopole field per- 
mits gauge transformations (in particular, the change of 
the direction of the spin quantization axis n and, hence, 
the positions of singularities). These do not change the 
equations of motion, but make a contribution to the La- 
grangian in the form of a total time derivative of some 
function of spin. This may in principle be significant for 
the calculation of the DW momentum. Naturally, the 
classical equations are not affected by the gauge trans- 
form. On the other hand, the term with dtp/dt in Eq. 
is of crucial importance for the quantum dynamics of do- 
main walls in spin chains with half-odd integer spins. 12 
For example, this term is responsible for the destructive 
interference of paths for DW tunneling from one mini- 
mum of the crystal potential to an adjacent one^ 

In order to construct a solution corresponding to a 
DW moving with velocity V, we have to write down 
the Landau-Lifshitz equations in angular variables, and 
to restrict ourselves to traveling wave- type solutions, 
9 = 0(£), <p = <p(Q, £, — x — Vt, with the natural bound- 
ary conditions, S x — > 5sin0o at £ — > ±oo, and 



S z -> aScos9 , S z 



-aS cos 8q 



at £ — > oo or £ — > — oo, respectively. Here a = ±1 is ttq, 
the topological charge of the DW and cos#o character- 
izes the ground state, see Eq. (J2J. Then, using Eq. J7J, 



the equations for 89/ dt = —V9' acquire the simple form 
v 9' = —2hsiiiip. Here and below the prime denotes dif- 
ferentiation with respect to £, and we have introduced 
the dimensionless variable 

v = V/V , V = aJS/h 

and put a = 1 in all intermediate equations. We will 
restore the dimensional parameters in final results only. 

The set of Eqs. ijHJ with the traveling wave ansatz has 
one integral of motion, which after taking boundary con- 
ditions into account can be written as 

i[(cos6>)'] 2 - sin 2 9 + 2hsin8 cos ip = h 2 = const. (12) 

Using the simple relation between 9' and sin ip, this can 
be rewritten in the variable 9 only. Finally we arrive at 
the simple differential equation for 9, 



(cos 9)' = crj \J2 sin 2 9 - v 2 + a 2 ^2h 2 -v 2 . 



(13) 



a i. 2 = ±1 are two independent discrete parameters 
which determine the topological charge of the kink as 
7To-topological soliton and the spin direction in the kink 
center, £ = 0, with respect to the magnetic field. Thus 
(Ti ! 2 fixes the type of the DW (i.e. DW with lower energy 
and DW with upper energy, see next paragraph). 

Eq. 1|13|) determines, in particular, the maximum pos- 
sible value of the velocity of the DW (7r-kink) , the critical 
velocity v c 



V2h 



(14) 



Eq. Ijl3(l can be integrated in terms of elementary 
functions, generalizing the solution for v = as given 
before i""' 2 ^ The analysis shows that the equation has 
both LDW and UDW solutions as described above. For 
v =/= the solutions are obtained by substituting 



cos 9 



i ip ■ yl — v 2 /2, sin a 



(7 2 



and can be given in the implicit form 



\f2o\ (£ — £o) = ip + tan a 1 ha 



sin[(V - a) /2] 



cos[(-0 + a) /2] 



(15) 



(16) 



As for v = Oii and in related models (see above) 
there are two types of DWs, with different energies. For 
v = 0, the DW with lower energy (LDW) in its center 
has the direction of spin S(0) parallel to the magnetic 
field H, whereas for the DW with higher energy (UDW) 
S(Q) is antiparallel to H. This applies similarly for 
moving DWs, however, 5(0) and H are not exactly 
parallel resp. antiparallel for V ^ 0. 

As for v = 0, the solution for the upper DW has a 
discontinuity in the space derivative 9'. At the critical 
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FIG. 3: The dependence of the DW energy (in units JS 2 ) 
and its momentum (in units TiS/a) on its velocity for the value 
of h = 0.5. 



velocity v = v c LDW and UDW become identical and 
the solution #(£) can be given in explicit form: 



cos# = cjcos#o ' sin(v2£/a) at |£| < a7r/2v / 2, 



and 



cos t 



a cos 6>o • (?/|£|) otherwise. 



It has a discontinuity in the second space derivatives. 

In order to find energy and momentum, the DW pa- 
rameters of interest, we do not need the solution for LDW 
and UDW in explicit form. It is easy to use Eq. (|13(l and 
to pass from the integrals over £ to the integration over 
9. Then an elementary calculation of the kink energies 
gives 



E, 



lower 



= JS 2 V2 



1- V - 
2 



arcsin , 



1-h 2 
l-v 2 /2 



En 



JS 2 ^{l-h 2 )(2h 2 -v 2 ) , 
Ei ower +4hJS 2 J(2-v 2 ) 



(17) 
(18) 



for the lower and upper DWs resp.. Thus, the depen- 
dence E(v) consists of two branches, Ei ower and E upper 
which merge at v = v c , and the full dependence E = E(v) 
is a continuous double- valued function, see Fig.|3 For the 
most interesting case of small magnetic fields, h <C 1, or 
v < v c <C vq, it is given by the unified equation 



E(V) = V2JS 2 



arcsin 



in v 7 ! - h 2 T 2/iWl 



(19) 



where the sign =p corresponds to the lower and upper 
branches of E = E{P) respectively. 

In the three-dimensional ferromagnet the UDW is un- 
stable. However, this instability develops as an inho- 
mogeneous perturbation (in the plane of the DW) and is 



not present for DWs in one-dimensional magnets. Below 
we will show that for the more natural representation of 
the DW energy, namely as a function of its momentum, 
the function E(P) is single- valued, and the upper branch 
just corresponds to the larger values of momentum. This 
explains its stability in the one-dimensional case. 

The kink momentum is determined as the total field 
momentum of the magnetization field, 23 , i.e. 



P 



d t s 4 s> 

5S 



-AdS. 

a 



(20) 



The dynamical part of the Lagrangian in Eq. @ and 
the expression l|2U|) for the momentum display singular- 
ities associated with the singular behavior of the vector 
potential A. The vector potential for a monopole in- 
evitably has a singularity on a line (Dirac string) and 
moreover is not invariant with respect to gauge trans- 
formations. The kink momentum also seems not to be 
invariant under these gauge transformations. If one uses 
a Lagrangian written in angular variables, problems due 
to the non-differentiability of the azimuthal angle ip at 
the points 9 = and 9 = 7r appear. (This problem is 
also important for the theory of moving two-dimensional 
topological solitons, see Ref. 28). It is important, how- 
ever, that the difference of the momenta for two differ- 
ent states of the DW is a gauge-invariant quantity^ To 
show this, let us imagine the DW as a trajectory in spin 
space {S}, S = S(£). The trajectories emerging from 
one point, say, S = S(e x s'm9 — e z cos9 ) and ending at 
another point S = S(e x sin9 + e z cos0 o ) can be associ- 
ated with DWs that move with different velocities but 
obey identical boundary conditions at infinity. In this 
case, the kink momenta are determined by integrals of the 
form J AdS along these trajectories. It is clear that the 
difference of the momenta is determined by the integral 
§ AdS along a closed contour. According to Stokes inte- 
gral theorem, the integral in question can be represented 
as the flux of the vector B = curlA through the surface 
bounded by this contour, J B dVL. Here the integral is 
taken over that region on a sphere which is bounded by 
the trajectories corresponding to the two kinks in ques- 
tion. The vector B = curL4 = S/\S\ 3 involves no sin- 
gularities. Returning to angular variables, this difference 
can be written as (US /a) ■ J sm9 dO dip, that is just the 
area on the sphere. Thus, the dependence V(P) or E(P) 
has been reconstructed, apart from the arbitrariness to 
choose the reference point for the momenta. 

The trajectories describing v = kinks appear to be- 
long to the large circle passing through the poles 9 = 0, ir 
and the end point of the magnetic field vector H . Ob- 
viously, the difference of the momenta P for the immo- 
bile LDW and UDW is determined by the area of the 
half-sphere, AP = 2nfiS/a. All the remaining trajecto- 
ries that correspond to moving kinks cover the regions 
between these trajectories. In particular, kinks at the 
critical velocities ±u c correspond to two non-planar tra- 
jectories. The v = ±v c trajectories can be reached by 
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proceeding from either type of v — domain walls, LDW 
or UDW. Thus, we can see that E(v) has two branches, 
one characterized by a higher energy and the other char- 
acterized by a lower energy. These two branches merge 
at the critical velocity v — v c . If we assume P = 
for the v = LDW, then the momentum grows up to 
P = ±P C , P c < 2irTiS/a as the absolute value of the kink 
velocity increases to v c . As we proceed further along the 
upper branch of the dependence E(V), the kink veloc- 
ity decreases, while the momentum increases further to 
P = ±2nhS/a when the velocity approaches again zero. 
When we continue, we begin to cover the area of the 
sphere once more and the momentum P grows, while the 
kink energy takes the same values as before: thus we do 
indeed arrive at a periodic dependence E(P) with period 
Pq determined by the total area of the sphere. The value 
of Pq depends only on the spin value S and the lattice 
spacing a, 



Po 



4nhS 



(21) 



The presence of two planar solutions (with ip = const) 
describing the DW's with v = and different ener- 
gies is a common feature for any model with uniaxial 
anisotropy subject to a transverse magnetic field. As we 
know, for Heisenberg exchange interaction and single-ion 
anisotropy, only numerical solutions describing mobile 
DW's can be founds But even in this case the peri- 
odicity continues to be the same. On the other hand, 
for the Ising model considered here, the exact solution is 
known, and the momentum can be calculated explicitly. 
As a result the momentum for LDW and UDW takes the 
following form 



2hS ( v I 1-h 2 

arcsin . ,. 

II \ V r : V 1 — IT / 2 



vhsVz 



h 2 



l~v 2 /2 



(22) 



Pupper — Plower 



4hS 
a 



1 • arcsin . 

2 



V 

- 23 



These equations together with Eqs. (|17l I18fl give us the 
dispersion law for the DW in implicit form as shown in 
Fig.H 

For the mostly interesting case of small magnetic fields 
h <C 1, or v < v c <C Vo the dependence P(v) can be given 
as the unified expression 



Pi 




± arcsin < 



(24) 



where ± correspond to UDW and LDW, respectively. Fi- 
nally, limiting ourselves to the case the dispersion 
law for both DW's takes the form 



E(P) = V2JS 2 



2h cos 

2 



2ttP 

~pT 



(25) 




FIG. 4: The dependence of the DW energy on its momen- 
tum for different values of the magnetic field (given next to 
the curves); arrows indicate the position of the momentum 
corresponding to the critical velocity v c . 



This equation will be used in the remaining part of 
the paper to describe the classical dynamics of domain 
walls and to perform their semiclassical quantization. It 
is based on the continuum approximation and this is the 
only conceivable approach to analytically describe mov- 
ing DW's; however, the validity of the continuum ap- 
proximation has to be justified before doing so since it 
seems to be non-adequate for the Ising model ferromag- 
net. This applies in particular to the limit h <C 1, when 
the DW thickness becomes comparable or even less than 
the lattice spacing a. To check the applicability of the 
continuum approximation we present in Fig. [5] a compar- 
ison of the DW energies as obtained from the numeri- 
cal approach to those obtained from Eq. (|25|l at P = 0. 
The discrepancy is seen to be of the order of 10% for 
h — ► and thus surprisingly small. In addition, the 
linear dependence of the DW energy on the magnetic 
field, E{H) = JS 2 (A±hB) is valid for both approaches, 
with values A w 2.0 and B w 2.0 for the numerical 
data and A = tt/V2 « 2.2 and B = 2^2 w 2.8 from 
Eq. (|25|l . As we will see below, the use of the dependence 
E{P) = JS 2 [A - hB cos(2n P/P )}, with "improved" co- 
efficients A and B, gives much better agreement with 
the numerical data for the dispersion law of quantum 
solitonsii^ 

One more important parameter, the maximum value 
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the parabolic approximation can be used and the energy 
takes the form E = P 2 /2M, where the effective mass of 
a kink is 



0.0 0.2 0.4 0.6 0.8 1.0 h 



FIG. 5: The dependence of the DW energy with zero velocity 
(in units JS 2 ) on the dimensionless magnetic field h; symbols 
depict the numerical data for the discrete model; the solid 
line is the theoretical result in the continuum approximation 
Eq. I125II : the approximation of numerical data by the linear 
function E = 2(1 — h)JS 2 is shown as dashed line. 
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FIG. 6: The ratio of the maximum value of the pinning 
potential to the quantity To characterizing the dispersion of 
the DW. For To, the "improved" value of the coefficient B — 2 
is used. Also shown as dotted line is the limit of this quantity 
for extremely small values of the magnetic fields—, Uo/Tq = 
0.5[1 -3(ft/2) 1/3 ]. 



M 



a 2 hj 



(26) 



It is seen from Eq. I|26[) that the effective mass turns to 
infinity at h — * 0. This is one more manifestation that 
in a purely uniaxial model of a ferromagnet in the ab- 
sence of a transverse magnetic field domain wall motion 
is impossible— I The use of the parabolic approximation 
(l2*f)|l seems to be adequate for DW's moving with small 
velocities. However, some important features are lost in 
this approximation: For example, the correct value of the 
energy bands forming the DW spectrum should be found 
within the analysis of the full periodic dependence E(P). 

We note that for small h the value of the effective 
mass for the upper DW branch, defined by E upper (P) = 
E(P /2) - {P-P Q ) 2 /2M up agrees with M from Eq. lf2%j) . 
However, M and M up are substantially different at finite 
h, see Fig. 0] above. 

To conclude this section, we emphasize that within 
the macroscopic classical approach the dispersion law 
(|25[l exhibits a periodic dependence of E(P) which, from 
Bloch's theorem, is characteristic for discrete quantum 
models. Moreover the maximum value of momentum, 
Pq = AirhS/a, which may be called the size of magnetic 
Brillouin zone, and the size on the crystalline Brillouin 
zone Pb — 2Trh/a have values of the same order at S ~ 1. 
The expression for Pq contains Planck's constant h and 
the lattice constant a and formally seems to character- 
ize both the quantum nature and the discreteness of the 
model. As discussed above, this is not the case: if one 
accounts for the macroscopic (continuous) character of 
the magnetization Mo = g^sS/a (otherwise the discus- 
sion of ferromagnetism is meaningless), Pq can be writ- 
ten using only Mq and the classical gyromagnetic ratio 
7 = e/2mc, Pq — AnMo/gj. 



IV. FORCED MOTION AND DISPERSION 
RELATION OF DOMAIN WALLS 



of the pinning potential Uq was treated numerically for 
arbitrary values of h. For extremely small values of h, 
U = E CB - E cs = 2JS 2 h[l - 3(/i/2) 1/3 ] was obtained 
from the analytical exp ressions for Ecs an d Ecb (see 
Eqs. (11,12) in Ref. Il4l). As we will see below, the ratio 
of Uq to the difference of the DW energies with P = 
and P = Pq/2, Tq = JS 2 hB is an important parameter 
for the description of the DW dynamics. Whereas the 
dependence of Uq on h is non-linear at small values of h 
because of the presence of the non-analytical term pro- 
portional to /i 1 / 3 , the DW "kinetic energy" To, is linear in 
h, see Sec.lVlbelow. The maximum value of Uq/T = 1/2 
is realized in the limit h — > 0, but in fact this important 
ratio decreases fast to rather small values when h grows 
(see Fig. |SJ). 

In the limit of small values of momentum, P — > 0, 



The periodic dependence discussed in the preceding 
section means that one can restrict the values of momen- 
tum P to the magnetic Brillouin zone, — Pq/2 < P < 
Pq/2. On the other hand, as with Bloch electrons, in 
order to analyze the motion under the influence of an 
external force, it is useful to consider that the DW mo- 
mentum obeys the equation dP/dt = F e , where F e is ex- 
ternal force, and to allow that P increases without limits 
beyond the first Brillouin zone. The expression for the 
DW energy E{P) can be used to describe the DW dy- 
namics in the spin chain in terms of the DW coordinate 
X considered as a collective variable. Its dynamics is 
governed by the Hamiltonian 



ttP. 



H(P, X) = EqQi) + T Q sin 2 ( — ) + U(X), 

-TO 



(27) 
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where Eq(K) is the DW energy at P = 0. Eo(h) has no ef- 
fect on the dynamics of the DW coordinate X and will be 
omitted below. Here we have introduced the parameter 

T = iV2hJS 2 = 2V2gn a SH 

which describes the magnitude of free DW dispersion, 
and added the potential U(X), without specifying its 
physical origin. For the discrete spin chain, the potential 
is the pinning potential originating from the discreteness 
effects introduced above. 

It is easily seen that the relation X = V — dH/dP 
as in analytical mechanics immediately gives the peri- 
odic (oscillating) dependence of the DW velocity on DW 
momentum 



vrTo . 2nP 
V = — sm (—) 

^0 P) 



(28) 



Inverting this equation we recover the expression for the 
momentum (1241 for small magnetic field. 

The dynamical equation for the Hamiltonian reads 
dP/dt = —dH/dX. Choosing different forms of the po- 
tential U — U(X) one can consider different problems 
such as the interaction between a DW and the inhomo- 
geneities in the medium or with external magnetic field 
H z directed along the easy axis. In the last case we have 



U{X) = -2 gi i B SH z ■ X/a, and P = 2g/j B H z S/a. (29) 

Therefore the DW velocity under the action of con- 
stant magnetic field (constant force) oscillates with time. 
Note that for a DW in the uniaxial ferromagnet this 
equation is nothing but the one-dimensional version of 
familiar Slonczewsky equations^! It describes the non- 
trivial properties of DW dynamics, such as the oscillat- 
ing motion of DW's as response to a constant external 
force. These effects were observed in a number of experi- 
ments on DW dynamics in bubble-materials, see Ref. l3ll 
Formally, such a motion corresponds to the Bloch oscil- 
lations well known for quantum mechanical electron in 
an ideal crystal lattice. Bloch oscillations for solitons 
in different media were recently reviewed by Kosevich. 32 
Also, such effects for DW's in spin chains with S =1/2 
have been discussed from the viewpoint of Bloch parti- 
cles recently^ However, it is important to note that for 
DW's in a continuum model the origin of the effects is dif- 
ferent in principle: It is not associated with discreteness 
and it exists even in the continuum limit. Thus, even in 
the classical continuum model the dynamics of a DW ex- 
hibits a number of properties peculiar to Bloch particles 
(electrons), i.e. to quantum objects moving in the peri- 
odic potential of the crystal lattice. On the other hand, 
the time dependence of this forced DW velocity is oscil- 
lating with the classical Larmor frequency £Il = 91 H z . 
Combining Eqs. I|28() and (|29fl one can find 

V(t) = V rnax sin Q L t, V m ax = ttTq/Pq. 



Thus, the quantum and classical regularities intertwine 
in a very intricate manner in the problem of domain wall 
dynamics. 

To analyze the quantum DW motion in the spin chain 
including discreteness effects, we will use Eq. 1)270 as a 
quantum Hamiltonian, using the periodic potential of 
Eq. ||IJ and Uq as found in the Sec. [H] Let us start 
with the case of nearly free motion, when Uq -C Tq. In 
this case the DW energy is given not by the usual mo- 
mentum, but instead by the quasi-momentum P which 
is determined only up to a reciprocal lattice vector. Ac- 
cording to Bloch's theorem, the energy E(P) should be 
periodic with period Pg — 2nh/a. When Uq — ► 0, the 
dispersion law is described by Eq. 1)25) 1. or Eq. 1)27) 1 with 
Uq = 0, for almost all P values (with small corrections 
~ U§). Only when the non-perturbed dispersion laws in- 
tersect, E(P) = E(P + nPs) with n integer, corrections 
become essential (~ Uq) and energy gaps appear. 

We now want to discuss the situation for different val- 
ues of spin S, since the ratio Pq/Pb = 2S depends on 
the spin value. Strictly speaking, small spin values, e.g. 
S = 1/2, 1, cannot be described in the frame of our 
semiclassical approach. But we argue that our approach 
will be valid at least qualitatively also for these cases. 

For spin S = 1/2 we have Po = -Pb, and there is no 
intersection of the non-perturbed dispersion law with its 
ima ge s hifted by Pb (the intersection of two parabolas in 
Ref. fl2l is an artefact of the parabolic approximation). In 
this case the effect of small U(X) on the DW dispersion 
relation E(P) is only quadratic in the small parameter 
Uq and is negligibly small. For all other spin values, 
S > 1/2, we have Po > Pb and dispersion curves Eq(P) 
extended periodically with period Pg = 2nh/a do inter- 
sect. However, all intersections occur on the boundaries 
of the Brillouin zone, namely, at P — ±Trh/a. It is clear 
that in this case the DW spectrum contains Nb = 2S 
energy bands. 



V. CLASSICAL DYNAMICS OF DOMAIN 
WALLS IN A FINITE POTENTIAL AND ITS 
QUANTIZATION 

In this section we will analyze the classical DW dynam- 
ics and its quantization for an arbitrary (finite) periodic 
potential of the form given in Eq. J3J. The most impor- 
tant ingredient to this analysis is that the Hamiltonian of 
Eq. (|2*7jl . just as the classical energy, is a periodic func- 
tion of momentum, and that there is an upper bound for 
the energy. As already mentioned, this property is al- 
ready present in the classical theory of DW motion and 
has nothing to do with quantum mechanics. We will show 
below that this property leads to results which are qual- 
itatively different from those obtained in the standard 
quadratic approximation, both for the pure classical DW 
motion as well as for its quantization. 

Let us consider the dynamical system described by the 
classical Hamiltonian of Eq. 1)27)1 . taking U(X) into ac- 
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count. The corresponding Hamilton equations 



dP 
~dt 



m dx 

dX' ~dt 



have an obvious integral of motion, 



nP 



&H 
dP' 



H(P, X) = T sin^ + U sin^ = const. (30) 

P a 

Here we have omitted the constant part Eo(h) which does 
not affect the equations of motion. These equations can- 
not be integrated in terms of elementary functions, how- 
ever, a sufficiently complete understanding of the DW 
dynamics can be found using the phase plane analysis. 
This system is characterized by the periodicity in mo- 
mentum P and by the presence of an upper bound for 
the Hamiltonian. In view of this its dynamics shows char- 
acteristic features which do not manifest themselves for 
standard dynamic systems with a parabolic dependence 
on momentum. 

It is easy to show that at arbitrary Uq/Tq there are 
two sets of center-like singular points in the phase plane, 
see Fig. One of them (CI— type center) corresponds 
to the minimum of the potential U(X) and the minimum 
of the "kinetic energy" T(P) = T sin 2 (nP/P ), 



PIP, 



1.0 



0.5 
0.0 
-0.5 

(a) 

1.0 
0.5 
0.0 
-0.5 

(b) 




PIP, 



1.0 XI a 




1.0 XI a 



CI : X = an, P = P m, H C \ = 0, 

the other set of centers (C2— type center) is located at 
the maximum of both potential and kinetic energies, 

C2 : X = a(2n+l)/2, P = P (2m+l)/2, H C 2 = U +T . 

Here and below n, m are integers. The existence of 
C2— type centers describing the steady small oscillations 
of the UDW near the potential maximum, is a unique 
property of Hamiltonian systems with an upper bound 
in the Hamilton function Ti(P,X). 

There are also two sets of saddle points in the phase 
plane, SI and S2. Points SI correspond to the maxima 
of the potential U(X) and the minima of T(P) 

SI; X = a{2n + l)/2, P = P m, H S i = U , 

points S2 correspond to the minima of U(X) and the 
maxima of T(P), 

S2; X = an, P = P (2m + l)/2, H S2 = T . 

At the allowed values of the integral of motion, Tici < 
H < Hc'2, the phase plane (P,X) separates into regions 
with different types of motion. Two types of finite motion 
are present within one interatomic distance (one period 
of the potential). One type of finite motion corresponds 
to the oscillations of the LDW near the minimum of the 
potential, it requires Hci < W < Hsi and is standard 
in the analytical dynamics of a particle. (We suggested 
Co < To here). A second type of finite motion, namely 
oscillations of the UDW near the potential maximum, is 



FIG. 7: The phase plane (P, X) for the DW dynamics for 
(a) Co = To and (b) Co — 0.25Tb- Separatrix trajectories are 
shown by full lines, other trajectories by dashed lines; full and 
open circles denote the positions of center-like singular points 
CI and C2, crosses correspond to the saddle points. 

realized for Hs2 < "H < Hc2- These regimes are sep- 
arated from the rest of the phase plane with different 
types of infinite motion by separatrix trajectories, end- 
ing in one of the saddle points. The only exception is the 
case Uq = To or Hsi — H-S2 when the separatrix trajec- 
tories connect saddle points of different type and form a 
two-dimensional net; as a consequence infinite motion is 
absent. For the transverse field Ising model with small 
Uo this is hard to realize, see Sec. [H] (in particular Fig.^JJ, 
but for methodological purposes we find it convenient to 
start with this case. 

In this case we have To = Uo, and one easily veri- 
fies that the area of the phase plane (mechanical action) 
per unit cell and one period of momentum is equal to 
Ao = 4arhS. Thus, from the Bohr-Sommerfeld quantiza- 
tion rules, there are (Ao/2Trh) — 2S quantum eigenstates 
of a DW, related to its oscillations per one unit cell. It is 
clear that if one takes into account tunneling transitions 
between equivalent points of the lattice, these localized 
states will turn into energy bands. 

This result coincides with the result obtained in the 
previous section for the opposite limit of an infinitely 
weak pinning potential. An additional argument on the 
number of energy bands is obtained from the well-known 
exact result for so-called Harper equation^ familiar in 
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the problem of electronic quantum motion in a periodic 
potential in the presence of a finite magnetic field. In ac- 
cordance with Ref. 33 the problem is reduced to a Hamil- 
tonian as in Eq. ^Dl, TLh = sin 2 (p/2) + sin 2 (7r/3g). For 
the rational case of the Harper equation, at /3 = m/n, 
the eigenvalue spectrum for the Hamiltonian Hh shows 
n non-overlapping bands. The simple canonical transfor- 
mation p — > irP/Po, q — > XPq/tt, leads to our Hamil- 
tonian H30fl with (3 — 1/2S. This gives immediately the 
above result, N b = 25. 

For T < Uq, a case more realistic for the Ising model, 
the situation is different: in addition to the localized tra- 
jectories in phase space, trajectories corresponding to the 
infinite motion above the potential barrier appear, see 
Fig. 03. These trajectories describe the overbarrier dy- 
namics of DW's in the pinning potential. It is clear, that 
such states should be described well by the model of the 
nearly-free particle, discussed in the previous section. 

The other limiting case, large Uo, i.e. Uo > To is hard 
to realize in the transverse field Ising model. On the 
other hand, it could be interesting for other models of 
ferromagnets supporting DW states, and we want to dis- 
cuss it briefly. The phase plane for large Uq/Tq can be 
obtained from Fig.[7|3 by replacing P/Pq by x/a and vise 
versa. Then, the topology of the phase plane is funda- 
mentally changed: trajectories with finite changes of the 
kink coordinate and infinite range in momentum, appear. 
For these trajectories the DW coordinate oscillates near 
certain positions, which do not coincide with a minimum 
or maximum of the pinning potential. This is nothing 
but Bloch oscillations in the pinning potential U(X). 

At this point we want to emphasize that the semiclas- 
sical result Aq = 4nhS does not depend on the details 
of the model. Only the periodicities of the Hamiltonian 
in P and X with periods Po and a are essential for the 
argument. The considerations presented in the previous 
section are model independent as well, and lead to the 
same value of Po for all transverse field models. Therefore 
the final result for the number of energy bands remains 
valid for more general models of ferromagnets subject to 
a transverse magnetic field. Thus the number Nf, of en- 
ergy bands for DW's in all these models should be equal 
to 25. This is in agreement with the numerical result for 
the kink dispersion law in the quantum transverse field 
Ising model with spin 5 = 5, where 10 energy bands were 
obtained^ 



VI. DW TUNNELING 

In the classical case, there are states corresponding to 
finite motion (oscillations) of DW's near the extrema of 
the potential in the phase plane both at small and large 
values of Uq. This applies to the LDW near the mini- 
mum of U(X) and as well to the UDW near the maxi- 
mum of U(X). Owing to the symmetry present in the 
Hamiltonian gJJ, P -> P /2 - P, X -> a/2 - X, their 
dynamics is described similarly; therefore one can con- 



sider only, say, the LDW case. In the quantum case, the 
tunneling of DW's from one site to another becomes pos- 
sible, and these states form the energy bands. It is clear 
that for extremely small Uo the nearly free approxima- 
tion discussed in section IIVI is valid. But the standard 
situation for semiclassical systems (like DW's for large 
spin ferromagnets) is the tight-binding limit, for which 
the probability of tunneling is small, and the states with 
energy E n are almost localized in a potential well. In 
order to estimate E n , we consider the parabolic approx- 
imation for both U(X) and T(P). In this case we have 
E n = Hloq ■ n, where 



h_ (<pu_\ yrmi 

W ° V M \dx 2 ) X=Q AhS ' 

M is the effective mass of Eq. <|26|) . Both Uo and To are 
proportional to 5 2 , therefore u>o oc 5 and for semiclassical 
spins S > 1 the value of E n for n <C 5 can be smaller 
than Uq even in the case Uq <C T). Then the width of 
the n th energy band, AE n , resulting from the tunneling 
between the quantized energy levels E n , is smaller that 
the value of To and even E n . It is clear that this first of 
all should correspond to the lowest level E n , i.e. n = 
(the "ground state of the DW" in the pinning potential), 
but it could be true as well for some higher levels with 
n > 1. 

In order to understand the possibility of semiclassical 
dynamics of the lower DW, consider the limit of small 
Uq, {Uq -C T)). The area Al under the separatrix tra- 
jectories 51 is easily found and the number of quantum 
states for the finite motion of the LDW, Nl = Al /2irh 
is defined by the following expression 



Nr. 



This is of the same order of magnitude as the value ob- 
tained from Nl ■ Tiluq = Uo. The same expression is ob- 
tained for the number of states of the upper DW, lo- 
calized near the potential maximum, Njj — Nl- At 
Uq <C To, the values of Njj and Nl are much smaller 
than the total number of DW states 25, but can still be 
large compared to 1: When the inequality 

-= < — <1, 

(meaningful for S> 1) holds, the states of finite motion 
are a small fraction of all DW states, but Nu — Nl ^> 1. 
It is easy to estimate the mean values of P 2 and X 2 in 
these states, 

(P 2 ) = P 2 

In view of that, it is not easy to satisfy the condi- 
tion of applicability of the parabolic approximation for 
both U{X) and T(P), (P 2 ) < (P 2 ) and (X 2 ) < a 2 
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simultaneously. Thus, in order to investigate the tun- 
neling splitting of the levels it is preferable to use the 
full non-parabolic Hamiltonian Eq. H30(l . The best way 
to analyze such models is the instanton technique, which 
can be applied for arbitrary field theories, see Ref. l34l 
In this approach, the amplitude of probability V12 of the 
transition of the DW from one given state to another is 
determined by the path integral J DX ■ exp{iA[X]/Ti}, 
where A[X] — J dt-C[X(t)] is the mechanical action func- 
tional. Here C(X, X) is the Lagrangian describing the 
dynamics of the DW coordinate X, and the integration 
DX goes over all configurations X(t) satisfying the given 
boundary conditions. The Lagrangian C{X 1 X) is easily 
obtained, as in standard mechanics, from C — PX — TL: 
express P and TL in terms of X, using the known expres- 
sions for P = P(V) and TL = E(V) + U(X) (see section 
llllll . replacing V by X. 

It is convenient to make a Wick rotation t — > ir, pass- 
ing to the Euclidean space-time. Then one has the Eu- 
clidean propagator J DX exp[— Ae/K\, where the Eu- 
clidean action is written as 

a e = J c E (x,n)dT, 

with Q = dX/dr. For the LDW this gives 



2HS 



Ce = — <>arcsinri ( ^-J + V c - y/V 2 + fl 2 +U(X). 

(31) 

Constructing the instanton solution can then be viewed 
as minimization of Ae with respect to X(t). The mini- 
mum of the Euclidean action is realized on the solution 
of the Euler-Lagrange equation for Ce- The first integral 
of this equation has the form 



2hS 
a 



(yv? + w-v c )-u{x) = c, 



with C = on the instanton solution. Then, the expres- 
sion for Ae reduces to the integral 

2hs r . , ,n,,„ 

Ae = / arcsmhf — )dX, 

a Jo V c 

which has to be calculated using the connection between 
n and X from the first integral, (QT /V C ) 2 = 4U{X)[T + 
U(X)]. Finally, the expression for the Euclidean action 
on the instanton solution takes the form 



A E 
% 



(32) 



VTo + U sin 2 V 
For the most natural case l7 <C T this gives the result 



A 



(o) _ 



8hS 



V at ^ « T o, 



(33) 



as usual for a Lagrangian with quadratic dependence on 
the velocity. For arbitrary Uq/Tq the calculations are 
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FIG. 8: The dependence of the Euclidean action on the value 
of Uo/To as found numerically (symbols) and from the asymp- 
totic Eq. J^J, A E oc ^Uo/Tq, valid for small U /T (dotted 
line). The maximum value of Ae/TiS = 1.533 is reached at 
Uo/To = 0.5, that is the maximum value of Uo/To for the 



transverse field Ising model. 



more complicated but the tendency is seen as follows: 
Ae is smaller than predicted by Eq. (f^Sf) for the same 
value of Uo/Tq, see Fig. |H1 Thus, the band width for 
finite Uq/Tq is larger than follows from Eq. (|33|l . 

For energy bands AE n with 1 « n « N L and U <C T 
we can use the fact that in the regime of classical motion 
the simple parabolic dispersion law is valid to calculate 
the tunneling splitting of the levels. This allows to use 
the standard expression for the tunnel splitting of the 
levels (see Ref. |3a, problem 3 after §50) and results in 



AE = (hujo l 1 ^) ex P ■ 



{2m[U(X) - E]} 1/2 dX/h 



Here a' is the turnover point of the classical trajectory 
with energy E n (corresponding to a non-split level) de- 
fined by the equation U(a') = E. It is convenient to 
rewrite it using the Euclidean action only. After some 
simple algebra we arrive at the following universal for- 
mula: 



AE„ 



hujo 



A E 

n 



-1/2 



exp 



A E 
h 



(34) 



where Ae is given by Eq. (|3"2"|l or Eq. (pT^ . with Ae oc 
S in both cases. This equation shows that the energy 
splitting for n > 1 continues to be exponentially small 
in S, but with larger pre-exponential factor. This is in 
good agreement with numerical dataiia 



VII. SUMMARY AND DISCUSSION. 

We start this summary with the explicit results ob- 
tained for the transverse field Ising model, and then go 
on to discuss the conclusions which might be obtained 
for more general models. The characteristic feature of the 
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Ising model is that the two important parameters Uq and 
T of the theory are governed by only one quantity, the di- 
mensionless magnetic field h; practically always Uq <C To 
holds. However, the behavior of the quantum DW spec- 
trum depends also on the value of the atomic spin 5. 
If Uo -C To, but S 2 > Tq/Uq, top and bottom energy 
bands are very narrow. The neighboring bands are nar- 
row also, but 5 timer wider, and the intermediate part of 
the spectrum consists of wide bands with narrow energy 
gaps. This picture corresponds well with the numerical 
data for the model with S = 5 of Ref. [IJ (see Fig. EJ. 

Unfortunately, the direct quantum mechanical ap- 
proach used in Ref. can be applied for small values 
of h only; then the discreteness effects are large and the 
continuum approximation is not appropriate. But the 
discrepancies are not too large, and even in this unfa- 
vorable situation we can claim at least qualitative agree- 
ment. The data are fitted quantitatively only with use 
of fitted dependence E{h) (see end of Sec. Inland Fig.|HJ). 
It is worth to mention that there are at least two rea- 
sons for discrepancies, caused by (i) nonapplicability of 
the continuum limit to the small field discrete classical 
Ising model, and (ii) the fact that the Landau-Lifshitz 
equation is not adequate to describe quantum DW's. We 
believe that discreteness is the main source of the dis- 
crepancy, and for models with wide enough DW's, (such 
as models of ferromagnets with nearly-isotropic exchange 
interaction and weak magnetic anisotropy), the descrip- 
tion of quantum dynamics of DW's based on the scheme 
proposed in this paper will be at least semi-quantitative. 

At this point it is interesting to compare the results 
of the instanton approach with data from an analysis of 
the model for extremely small field, when spin zero-point 
fluctuations far from a DW can be neglected and the one- 
site approximation can be appliedpM In both approaches 
the widths of the energy bands with lower energies are ex- 
ponentially small in the spin magnitude 5. The one-site 
approximation yields an exponential factor of the form 
of exp(— 251n2) w exp(— 1.395), 14 which is quite close 
to the value exp(— Ae/K) — > exp(— 1.535) at h — ► 0, ob- 
tained in the instanton approach using improved data 
for T , C/o/To — * 1/2. It is interesting to note that using 
the data obtained directly from the continuum approx- 
imation, Uq/Tq = l/(2v2), the factor takes the form 
cxp(— 1.345) which is much closer to the value obtained 
in Ref. The agreement between the data can be con- 
sidered rather good, regarding the conceptual difference 
between the two approaches. 

Good agreement is also found with respect to the fact 
that the band width increases with energy, the ratio of 
the widths of neighbouring energy bands is given by a 
factor of 5. But the pre-exponential factors for these, 
essentially alternative, approaches differ by a factor of 
V5. It is most likely that taking into account fluctua- 
tions far from the DW center becomes essential for the 
determination of the pre-exponential factor. It is worth 
noting the formal agreement that, within the instanton 
approach, this factor is determined by the fluctuation 




FIG. 9: The dependence of the DW energy on its momentum 
for the magnetic field h = 0.05 as obtained in the continuum 
approximation. The dashed line presents the result of the 
continuum approach, the solid line is improved by using the 
fitted dependence E(h) for small h (see Fig. Symbols 
show the numerical data for the dispersion law of the quantum 
DWpii 



determinant only, again, far from the instanton center. 
Again, the compliance of data of these two approaches 
can be considered as satisfactory. 

The instanton analysis demonstrates that while the 
field is increasing, Ae/Ti very rapidly decreases, see 
Fig- El Thus, the width of the energy band grows rapidly, 
the gaps become narrow and the nearly-free limit can be 
realized for all energies. At 5 — 5 the value of Ae/^ 
becomes of the order of unity and AEq ~ Tiujq already at 
Uq/Tq < 0.008 (equivalent to the magnetic field becom- 
ing close to h ~ 0.3). But for any values of the parameter 
h, 25 energy bands are present in the DW spectrum. As 
we argued above, this result is not specific to the Ising 
model. It is mostly due to the value of the period P , 
which is universal for all models of uniaxial ferromagnets 
in a transverse magnetic field. Thus, the result that 25 
energy bands exist should be valid for all transverse field 
models of ferromagnets. 

The results obtained here for ferromagnets can be 
applied to Ising antiferromagnets as well. The simple 
canonical transformation 



GUy 



(-1)' 



Oil 



results in the same commutation relations for staggered 
magnetization a n as for usual spins 5„ and reproduces 
the same form for the Hamiltonian, Eq. replacing 
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S n by a n and changing the sign of J. It is easy to 
prove, after the appropriate gauge transformation, that 
the form of the Lagrangian, Eq. H10[) is also unchanged. 
Thus, all the results obtained here for the DW dynamics, 
its dispersion law and its semiclassical quantization are 
valid for Ising antiferromagnets as well. Note that such 
a simple connection can be true only for models with a 
Hamiltonian depending on two components of spin only. 
This excludes near Heisenberg antiferromagnets. For the 
near Heisenberg antiferromagnets with weak anisotropy 
of any origin, exchange or single-ion, treated on the basis 
of the sigma-model, the singularities in the Lagrangian 
for the staggered magnetization are absent for a suffi- 
ciently general set of models (including antiferromagnets 
subject to strong external magnetic fields and in the pres- 
ence of Dzyaloshinski-Moriya interactions of arbitrary 
symmetry^) and a periodic DW dispersion law is not 
realized^ 

Another interesting feature, which is common to clas- 
sical and quantum physics, is that states with finite range 
of motion for the coordinate, but infinite range for mo- 
mentum are present for large enough pinning potential, 
Uq > Tq. Evidently such trajectories do not exist for 
standard mechanical systems with a parabolic dispersion 
law H = P 2 /2M + U{X). Of course the number of quan- 
tized states of DW motion per unit cell does not depend 
on the ratio of and To and remains equal to 2S. For 
the case of extremely weak dispersion T -C Uq, it follows 
that intermediate DW positions in the periodic potential 
are quantized and are equal to 2S discrete values. Thus, 
in the case of semiclassical quantization of DW motion 
within the continuum problem, there are 2S positions of 
the DW, i.e. some internal scale, a m — a/2S appears in 
the problem. This fact might be used to explain intu- 
itively that a magnetic Brillouin zone, Pq = 2SPb, ap- 
pears (here Pb = 2nTiS/a is the size of the Brillouin zone 
in a chain with interatomic distance a). To explain this 
in a slightly different way, we mention that the DW coor- 
dinate is directly connected with the z— projection of the 
total spin S z . tot, see pj, and this quantum-mechanical 
quantity can change by unit steps only. 

It is interesting to compare the features established 
here for transverse field models to the 180-degree DW 
(n— kink) in ferromagnets with rhombic anisotropy as 
considered by Braun and Lossi^. For their model there 
are two different DW's with lowest energy (Bloch DWs) 
and the corresponding trajectories are the halves of 
the large circle passing through the easy axis and the axis 
with the intermediate value of anisotropy. 

It is easy to argue, that in this case the dependence 
E(P) found in the continuum approximation is again 
periodic, but the period of this dependence Pq^ is de- 
termined by half of the sphere S 2 = S 2 , and it is two 
times smaller than that for the transverse field model, 
P a ,n = 27rn5/ai 17 i 23 i 29 The specific model treated iiM is 
exactly integrable, and the free dispersion law E(P) with 
this period can be constructed explicitly, see Refs. 00- 
On the other hand, for the full classification of DWs, i.e. 



to describe the presence of two different, but energeti- 
cally equivalent DWs one can in this situation introduce 
the discrete parameter chirality. Braun and Loss have 
shown from a different argument that the minima of the 
dispersion laws E(P) for DWs with different chiralities 
have to be located at different points of the P-axis, with 
distance equal to Pq^ (in our terminology). Owing to the 
parabolic approximation, used in Ref. Il2t the dispersion 
curves E{P) for different chiralities seem to be intersect- 
ing in some point. But accounting for the real periodic 
dependence for E{P), these two minima are smoothly 
connected, with the maximum at the energy of the Neel 
DW, and no intersection occurs. Thus, there is no reason 
for the chirality tunneling to take place in this formula- 
tion. The deep connection of chirality and linear momen- 
tum and the suppression of quantum tunneling of chiral- 
ity for free 180-degrees DWs in ferromagnets has been 
discussed recently by Shibata and Takagipi^ It is worth 
to note here, that also for 180-degree DWs in Heisenberg 
antiferromagnets two non-equivalent DWs with equal 
energies are present, but the dependence E{P) is not pe- 
riodic. Chirality is not connected with the value of linear 
momentum and the effects of quantum tunneling of chi- 
rality are present even without the pinning potential. 8 

Finally we want to discuss the perturbation of the dis- 
persion law for a 180-degree DW caused by a small pe- 
riodic pinning potential U(X). For this case, one has to 
locate the intersections of non-perturbed spectra E(P), 
periodically extended with period Pg = 2nh/a. The sit- 
uation here is slightly more complicated than for trans- 
verse field models: For small spins S = 1/2 and S = 1 
one has Pq^ < Pb, intersections are absent and no gaps 
are present (this is the same situation as for the trans- 
verse field model with spin S = 1/2). For integer spins 
S = k, k > 1 intersections occur and energy gaps appear 
at the intersection points. However, all these gaps are 
located on the boundaries of the Brillouin zone. Then, 
the number of energy bands is equal to S. But for suffi- 
ciently large half-integer spins S — (2k + 1) /2, k > 1, one 
has Pq — Pe(2n + l)/2, and some intersections are inside 
the Brillouin zone at P = (2n + l)nh/2a, with integer 
n. This leads to the appearance of some additional gaps, 
which could be associated with the tunneling of chiral- 
ity. In this case the number of gaps increases and for 
half-integer spin S the number of energy bands in the 
dispersion relation for 180-degree DWs is equal to 2S. 
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